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Abstract 

Type la snpernovae have played a crncial role in the discovery of the dark 
energy, via the measurement of their light curves and the determination of the 
peak brightness via fitting templates to the observed lightcurve shape. Two spec¬ 
troscopic indicators are also known to be well correlated with peak luminosity. 

Since the spectroscopic luminosity indicators are obtained directly from observed 
spectra, they will have different systematic errors than do measurements using 
photometry. Additionally, these spectroscopic indicators may be useful for stud¬ 
ies of effects of evolution or age of the SNe la progenitor population. We present 
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several new variants of snch spectroscopic indicators which are easy to antomate 
and which minimize the effects of noise. We show that these spectroscopic indi¬ 
cators can be measured by proposed JDEM missions such as SNAP and JEDI. 

Subject headings: cosmology: dark energy — stars: atmospheres — supernovae 


1. Introduction 

Type la supernovae are now recognized as one of the most important cosmological 
probes. Astronomical interest has been focused on Type la supernovae since it was recognized 
long ago (Wilson 1939; Kowal 1968) that they are good “standard candles” and hence are 
useful cosmological probes. Their use as cosmological beacons led to the discovery of the 
“dark energy” (Riess et al. 1998; Garnavich et al. 1998; Perlmutter et al. 1999) and they are 
an integral part of plans to further characterize the nature of the dark energy equation of 
state with a wide held satellite probe, the “Joint Dark Energy Mission” (JDEM) as well as 
ground-based studies. The reliability of SNe la as distance indicators improved signihcantly 
with the realization that the luminosity at peak was correlated with the width of the light 
curve (Phillips 1993) and hence that SNe la are correctable candles in much the same way 
that Cepheids are (Phillips et al. 1999; Goldhaber et al. 2001; Riess et al. 1995). To date, 
cosmology with Type la supernovaehas used photometry exclusively over spectroscopy, except 
that a spectrum is required for supernova type conhrmation and can be used for the redshift 
determination, the luminosity at peak is always determined photometrically (by matching 
lightcurve templates). However, there exist spectroscopic luminosity indicators hrst discussed 
by Nugent et al. (1995). Garnavich et al. (2004) showed that iHsi (Nugent et al. 1995) is well- 
correlated with Amis. Of course, the value of Mb at maximum light is also well correlated 
with Amis (Phillips et al. 1999; Phillips 1993). Thus, the spectroscopic luminosity indicator 
should be well correlated with the peak luminosity, making it a useful secondary luminosity 
indicator for cosmological studies. Even though it is correlated with peak luminosity via the 
photometric indicator Amis, it is still useful as an independent check because it will have 
different systematic errors than light curve shape htting via templates. 

We study the correlation of the two spectroscopic luminosity indicators ifisi and 
with luminosity, both for observed supernovaeand for synthetic spectra of the parameterized 
explosion model W7 (Nomoto et al. 1984) calculated using PHOEMIX(Hauschildt & Baron 
1999, 2004) in LTE. We dehne new, more robust spectroscopic luminosity indicators and we 
show that these indicators can be used with good accuracy via proposed JDEM missions 
such as SNAP (SNAP Gollaboration 2005) and JEDI (Crotts et al. 2005). 
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2. Spectral Indices ^si ^ ^Ca 

While most Type la supernovaefall into the category of “Branch normal” (Branch et al. 
1993) with B — V colors at peak close to 0, there are known to be a class of bright ‘T991T-like” 
supernovaethat have weak Si ll lines at maximum, but which develop characteristic Si ll lines 
shortly after maximum light. Additionally there is a class of dim “1991bg-like” supernovae, 
which are quite red and are thought to produce very low amounts of nickel, ~ 0.1 M© 
as opposed to ~ 0.6 Mq for Branch normals (Mazzali et al. 1997). 

Nugent et al. (1995) showed that the diversity in the peak luminosity from the dim 
1991bg-likes, to the bright 1991T-likes could be interpreted in terms of a single parameter 
^modeB^ synthetic spectral models, with the dim, low nickel mass 1991bg-likes being cool 
and the bright 1991T-likes being hot. Furthermore, Nugent et al. (1995) found that the ratio 
of the depth of the characteristic dehning feature of SNe la, Si llA6355, to that of another 
silicon line Si IIA5972 is correlated with the spectral sequence and hence with the peak 
luminosity. The obvious interpretation is that this is a sequence in nickel mass produced, 
although that interpretation is not completely proven. Nugent et al. (1995) also dehned 
another spectroscopic indicator using the ratio of two peaks near the Ca II H+K feature. 


2.1. Definition of ^si 

^Si is named for the characteristic Si ll trough found at ~ 6100 A and usually associated 
with the absorption trough of the Si llA6355 Si ll line P-Cygni feature. 

Fig. 1, illustrates the dehnition of 3^5^, using the observed spectrum of SN 1992A 3 days 
after maximum. The dashed lines are the polynomial hts used to determine the maxima 
found in the wavelength zones dehned in Table 1. The two depths are indicated by arrows. 


Table 1. dehnition zones 


Feature 

Lower Wavelength (A) 

Upper Wavelength (A) 

First maximum 

5500 

5700 

Si IIA5972 minimum 

Variable 

Variable 

Second maximum 

5850 

6000 

Si llA6355 minimum 

6050 

6250 

First maximum 

6200 

6450 



arbitrary units 
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Wavelength (A) 

Fig. 1.— The definition of the method used for determining ^si- The observed spectrnm is 
that of SN 1992A 3 days after maximum and the dashed lines indicate the polynomial fits 
to determine the maxima and minima. The arrows indicate the depths and used 

to compute 
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Fig. 2.— The definition of ^ca- The dashed lines indicate the polynomial fits used to 
determine the maxima, which are indicated by the arrows. 
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The Si llA6355 minimum band is centered on the Si llA6355 line blueshifted approx¬ 
imately 10,000 km s“^(Due to homologous expansion, it is customary to use velocities 
interchangeably with wavelength shifts when analyzing supernovaespectra.). The value of 
10,000 km s“^corresponds to the average velocity at which the Si ll line is observed to form. 
The ~ 200A wide wavelength zone where we search for the second minimum is noted as 
“variable” because it is centered on the Si llA5972 line blueshifted by nearly the same actual 
velocity as the Si llA6355 minima. This to ensure that the right trough is found, which is 
assumed to be at least partly created by the Si llA5972 Si ll line, thus it should form in the 
same physical region. 

is dehned as the ratio of the vertical distances between each minima and the reference 
lines drawn between two consecutive maxima and as shown in Fig. 1, 

Ss. = ^. (1) 

Ored 

This dehnition differs slightly from the one of Nugent et al. (1995), where tangents to 
the spectrum were manually selected instead of the reference lines between two consecutive 
maxima. These two methods mainly differ when the zones searched for maxima only display 
inflection points, but our method was easier to automate. 

Finally, because it is dehned as the distance between points of the spectrum and ref¬ 
erences linearly coupling points of the spectrum, 3?5iis independent of the absolute hux 
calibration. Moreover, as it is calculated within a wavelength region only ^ lOOOA wide, it 
will be relatively insensitive to both the relative hux calibration and reddening. 


2.2. Definition of 3?ca 

Following Nugent et al. (1995), we dehned two zones in the 3000A-4000A region of the 
spectrum, centered on the Ca ll 3650A and 3933A lines, as shown in the right hand side of 
Fig. 2. ^ca is dehned as the ratio of the ~ 3933A maximum, max(F(3933)), to the one at 
^ 3650A, max(F(3650)). 

^ max(F(3933)) 
max(F(3650)) 

The width of the zones searched for these maxima are displayed in Table 2. 


The maxima are almost never found at the exact wavelength of the Ca ll lines. This is 
expected since even if we had pure Ca ll, effects other than resonant-scattering could move 
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the feature’s maximum wavelength. Moreover, these features are line blends, and not pure 
Ca II lines. 


The study of these wavelength zones with PHOENIX showed that the dominant elements 
of these blends had lines displaying the same qualitative evolution as ^ca with luminosity. 
This leads us to propose that the maxima ratio correlation with luminosity is only a corollary 
of the existence of a correlation between the blends. Hence, we defined an integral ratio we 
named ^cas (the S stands for surface) as 





( 3 ) 


The wavelength regions were empirically chosen in order to group all the contributions with 
the same evolution with luminosity than 3?ca- We chose the wavelength limits in Eq. 3 
in order to produce a spectral index that is well correlated to luminosity and robust to 
wavelength errors. 

Both the ^ca and ^caS ratios, like ^si, are independent of the absolute flux calibration, 
and not too sensitive to the relative calibration quality or reddening. 


2.3. Dealing with Noise: 

The main merit of SfJc'aS over ^ca is the increased signal to noise ratio. On the other 

hand, 3?ca and SfJs'j will be highly sensitive to noise since they rely on maxima and minima 

fh 

detections. In order to reduce this effect, we fit each wavelength zone considered with a 4*^ 
degree polynomial and search for their local maximum or minimum. 

Because we want and ^ca calculation algorithms to be automatic and robust, we 
need to search wavelength regions wide enough to be sure to find the right maxima and 
minima. But in these large zones the spectral feature shapes are more complex than that of 
a second degree polynomial. On the other hand, noise may generate oscillations which will 


Table 2. ^ca Maximum Zones 


Central Wavelength (A) 

Zone half-width (A) 

3650 

48 

3933 

65 
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be followed by a polynomial of too high a degree. We chose 6^^ degree polynomials as they 
empirically proved to be a good compromise between these two effects. 


3. Correlation of 3?ca with Luminosity 

We applied these definitions to the public supernovaespectra we were able to gather, 
listed in Table 3. The absolute magnitudes of Reindl et al. (2005) are derived for a common 
arbitrary value of the Hubble Constant Hq = 60km s“^ Mpc“^. The values of the reddening 
and (5mi5 were chosen to minimize deviations from the Hubble law for a set of 66 SNe Iain the 
Hubble flow. While other prescriptions for determining Mb would give somewhat different 
results, with our small sample it is advantageous to use a single consistent prescription to 
determine Mb rather than relying on a variety of methods to determine the distances to 
nearby SNe la. When present surveys that will construct large homogeneous samples of 
SNe Iain the Hubble flow are complete, it will be simple and necessary to redo our analysis. 
Nevertheless, we don’t expect to alter our qualitative conclusions. Reindl et al. (2005) were 
unable to determine a unique solution for the reddening to SN 1986G, and thus we have two 
Mb values for SN 1986G which we carry through in our analysis. In the following, we will 
calculate results for both values of the absolute B magnitude of SN 1986G, we will refer to 
the Mb = —17.58 case as “Gase A”, and —19.16 as “Gase B”. Since the analysis of Reindl 
et al. (2005) is the most recent and careful, we consider Case B to be the more likely value 
for Mb and we will present plots of our results for Case B only, however since the reddening 
correction is uncertain and we are dealing with only a small number of supernovae, we will 
present our results for both cases in tabular form. For our calibration, we restrict ourselves 
to “Branch Normals” which we define to be SNe lawith —19.0 < Mb < —19.8, thus in Case 
A, SN 1986G is not included in the set of calibrators, whereas in Case B it is. As more 
spectra of SNe la with well-measured values of Mb become available they are easily added 
to our analysis. 

Even though iftsi and iftca were originally defined for SNe laspectra at maximum light, 
we used spectra ranging from 5 days prior to maximum to 5 days after maximum. We 
examined various epochs around maximum in order to determine whether or not the date of 
the spectra has a critical impact on the ratio’s correlation with luminosity. Moreover, since 
space programs like SNAP will only take one spectrum per supernova, we wanted to quantify 
the time window where spectra are usable for luminosity measures with line ratios. Note 
that JEDI proposes to take multiple spectra of each SN, at 5-day intervals in the observer’s 
frame. 

On the other hand, because of time dilatation, the rise time for supernovaeat ~ 1 is 
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twice as slow in the observer’s rest frame. Relaxing the time constraint to ±5 days aronnd 
maximum was thus conservative. 


3.1. 3?ca Correlation with Luminosity 

Some of the spectra listed in Table 3 lacked the wavelength coverage needed to calculate 
SfJc'a and ^cas, those epochs are indicated with an asterisk. Neglecting these, the results for 
Case B are plotted in Figs. 3 and 4. 

The 3?ca values for maximum light spectra agree with Nugent et al. (1995) to within a 
few percent. In order to calibrate the luminosity relations, we computed the linear regression 
of 3?ca and 3?cas values with respect to blue magnitude. This linear regression should be valid 
in the restricted set of “Branch Normals”. 3-D deflagration models have difficulty reproducing 
the very low nickel mass of SN 1991bg-like supernovaeand we illustrate below that synthetic 
spectra of W7 (Nomoto et al. 1984) indicate that there may be a departure of the spectral 
ratios from linearity for 1991T-like supernovae. Thus, by restricting our calibrators to Branch 
Normals we might be using a more similar population. From the linear regression we obtain 

^Ca = Mb + (4) 

^CaS = ^ ^^CaS 

We list the resulting parameters in Table 4. Since we have no information on the 
experimental errors, we set them arbitrarily equal to unity. 

We also calculated the standard deviation oso for ^ca and 3?ca5, defined as the 
quadratic mean of the distance of each supernova from the linear regression, and have listed 
them in the same table. This definition of the error combines measurement errors, time 
dependence, and intrinsic dispersion. 


Luminosity measure precision using ^ca or With the linear regression cal¬ 
ibration of ^ca and with respect to blue magnitude and the associated variance, we 

where able to estimate the luminosity precision measure using these ratios as follows: 


^Ca 

AL 


= 




Ca 


O'Mr. 


^Ca 


tCa 




Ca 




Ca 


L 
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Table 3. List of Observed supernovae 


SN 

Epoch (day wrt 

max) 

Ms 

SN 1998aq 

-3, 0, 1, 2, 3, 4 


-19.56 

SN 1981B 

0* 


-19.54 

SN 1998bu 

-4, -2, -1 


-19.49 

SN 1994D 

2% 3*, 4*, 5*, -3, 

-4, -5 

-19.47 

SN 1996X 

-4, -2, 4, 0% r 


-19.43 

SN 1989B 

-1,-5 


-19.42 

SN 1992A 

0, -1, -5, 3, 5 


-19.32 

SN 1986G 

-1, -3, -5% 1, 3 


-19.16 or -17.58 


’^Supernovaeused in our simulation. Absolute magni¬ 
tudes are taken from from Reindl et al. (2005). The 
spectra from the epochs denoted with an asterisk lacked 
coverage in the blue and were excluded from the calcu¬ 
lations of 3?ca and ^caS (see text). 


Table 4. ^ca & luminosity measure precision 






3ftca (Case B) 

0.20 

2.1 

0.10 

(Case B) 

0.31 

3.1 

0.10 

Sica (Case A) 

0.19 

1.2 

0.16 

^Cas (Case A) 

0.39 

1.4 

0.28 







RCa 
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Fig. 3.— 3fJcaCalculated using supernovselisted in Table 3 (Case B). The line is the result of 
the linear regression. 



RCaS 
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Mb 


Fig. 4.— 3fJcascalculated using supernovselisted in Table 3 (Case B). The line is the resnlt 
of the linear regression. 
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We display in Table 4 the present measnrement accnracy on the bine magnitnde as 
estimated for ^ca and 3?caS- We hnd in Case B that the accuracy of the method reaches 
the 0.1 magnitude level. 


3.1.1. ‘Stsi correlation with luminosity: 

The correlation calibration Fig. 5 shows the iksi values calculated with supernovaefrom 
Table 3, including spectra at all epochs. We excluded a negative value from the —5 day 
SN 1994D spectrum, for which the blue trough associated with the Si llA5972 feature does 
not exist. The SfJs'j ratio is of course mathematically well dehned even in this peculiar case, 
but since this point falls out of the general trend and increases the 3^5^ slope, we removed 
it from the linear regression. This selection is compatible with observational constraints, as 
such spectra are easily identihed when the signal to noise is large enough to measure 

The linear regression and the corresponding relative dispersion for this selection of 
supernovaeis also shown in Fig. 5. 


We also calculated the associated standard deviation and found it to be ~ 0.09 
(Case B) or ~ 0.06 (Case A). This variance includes the measurement errors, time 

dispersion, as well as the intrinsic dispersion of the supernovae. 


Luminosity measure precision: As in §3.1 we estimate the accuracy of the luminosity 
determination using directly from the dispersion with respect to the linear regression. 
Table 5 summarizes the results. 

Table 5. luminosity measure precision 






O'Mb 

3^5* (Case 

B) 

0.05 

0.75 

0.07 

ihtsi (Case 

A) 

0.06 

0.59 

0.10 
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- 19.0 


Fig. 5.— Sfts'jcalculated using supernovsein Table 3 (Case B). The line is the result of our 
linear regression. 



arbitrary units 
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Fig. 6.— The definition of the method used for determining ^siS- The arrows indicate the 
maxima used to compute ^siS- 
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4. A New Spectral Ratio 


Due to the correlation with luminosity, we have focused on the formation of features 
in the zone. More will be said about the line formation in this spectral region in a 
forthcoming publication (S. Bongard et ah, in preparation), but the clarihcation of the 
spectrum formation process led us to devise a new line ratio indicator which we call ^siS- 

^Sis is dehned as the ratio of two maxima. The hrst maximum is the one used in the 
calculation, located around Si llA6355, and the second one is the bluer peak of due to 
sulfur (hence the name Si+S a mixture of silicon and sulfur) located in the 5600 — 5700A 
region as shown in Figure 6. The zones where we look for these maxima are listed in Table 6. 
As was the case for ^si and 3?ca, we used a 4^^ order polynomial to ht the regions searched 
for the maxima, in order to increase robustness of the ratio with respect to Poisson noise. 


We also dehne the integral ratio ^siss as: 






5700 

5500 


/< 


6200 

6450 


FxdX 

FxdX' 


(7) 


Once again, the simplicity of this integral ratio calculation and its insensitivity to noise 
makes it an interesting alternative to the line htting method. 


5. correlation with luminosity 

We plot the linear regression of ^sis correlation with luminosity in Fig. 7 We also plot 
in Fig. 8 the ^siss correlation with luminosity. Table 7 summarizes the results for both 
quantities. 


Table 6. zones 


Lower Wavelength (A) Upper Wavelength (A) 

First maximum 6200 6450 

Third maximum 5500 5700 
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Fig. 7.— ^sis (Case B) calculated using the supernovseof Table 3. 
our linear regression. 


_I_ 

- 19.0 


The line is the result of 
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Fig. 8.— ^siss (Case B) calculated using the supernovseof Table 3. The line is the result of 
our linear regression. 
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5.1. Luminosity measure precision 

We summarize the different slopes and dispersion valnes for ^sis and ^siss calcnlated 
with or withont the time dispersion correction in Tables 7 and 10. 

With this spectral indicator family, even without time correction, we have augmented 
the blue magnitude accuracy by a factor of two compared to ^ca- In Case A, SfJs'j and 
have similar accnracies; however, in Case B ^sis is slightly better. Since ^sis and cannot 
be considered as independent, we still have only two independent spectral indicators with 
^Ca- But since one of them is as accurate as the light curve method, coupled with enough 
supernovaethey have the possibility to constrain evolutionary effects. 

6. Time Evolution of Spectral Indices 

We plot in Fig. 9 and Fig. 10 the ^ca and 3?ca5 time dependence for each snpernova. 
Whenever enongh data points were available, we interpolated the t = 0 point nsing a 2*^^ 
degree polynomial. For the cases with two points only, we used a straight line. 

The snpernova SN 1994D had enongh points for a 2*^^ order polynomial interpolation, 
bnt since all of them are for t < —3 days, and since the trend is different from the other snper- 
novae, we considered the t = 0 point to be too mnch of an extrapolation to be nsed. Whether 
the SN 1994D dispersion in ^ca or ^cas is intrinsic or due to measurement uncertainties is 
still an open question. 

Replacing each ^ca and ^cas by the corresponding t = 0 calcnlated valne, we recom- 
pnted the linear regression, and the associated standard deviation. These resnlts are shown 
in Table 8. 

More spectra are needed in order to show if there is a correlation in the evolntion ^ca 
Table 7. ^sis Inminosity measnre precision 



(7 

a 

(^Mb 

3^5*5 (Case B) 

0.06 

0.67 

0.09 

'^SiSS (Case B) 

0.06 

0.73 

0.09 

'Sisis (Case A) 

0.06 

0.43 

0.14 

^siss (Case A) 

0.07 

0.53 

0.11 



Fig. 9. 
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Fig. 10.— 3?caS time evolution 
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or with epoch, but it is apparent that there is no uniform trend with time, and the 

attempted correction to maximum light has increased the dispersion from Table 4 to Table 8. 
Therefore, it is conservative to treat the time dependence of 'Slca and 3?cas as an uncontrolled 
dispersion. 

Fig. 11 displays the 'Slsi time dependence of the supernovaewe used. We calculated a 
quadratic regression for each supernova with three or more values. We then replaced 
each value by the corresponding value at t = 0 obtained via either interpolation or 
extrapolation. The new linear regression calculated for 3^5* correlation with luminosity is 
also plotted in Fig. 11. We calculated the same correction for the less favored luminosity 
choice for SN 1986G (Case A). 

For the favored SN 1986G luminosity case (Case B) the extrapolation to t = 0 sig- 
nihcantly improves the measurement, as seen in Table 9 and lowers the dispersion to 0.1 
blue magnitude so that the ratio probes “Branch normal” SNe laluminosities with good 
accuracy. 

The standard deviation calculated with respect to this new regression, still removing 
the events with negative is now ~ 0.07 (Case B) or ~ 0.12 (Case A). 

The top panel of Fig. 11 shows that there is no general time dependence trend for our 
supernovaesample. 

In Fig. 12 and Fig. 13 we plot the time dependence of and 'Slsiss as well as the linear 
regression calculated on the values at t = 0 interpolated or extrapolated with quadratic hts. 
The three supernovaewith sufficient data do not display a common time evolution, making 
doubtful the existence of an universal ‘Stgis or ‘Stsiss time correction. 

We again calculated the standard deviation, the results of the time correction are sum¬ 
marized in Tab. 10. 

Correcting for the time evolution increases the '^siS and '^siss efficiency both by in¬ 
creasing the slope and decreasing the dispersion, the only drawback is the need of at least 
three points to interpolate or extrapolate to t = 0 with a quadratic ht. Since there is no 
common trend for the three supernovaethat were ht, this correction remains purely empirical. 


7. Comparison with synthetic spectra 

In order to ascertain whether the trends that we have found for the spectral indices 
can be reproduced by models, we have calculated a set of local thermodynamic equilibrium 
(LTE) models using the generalized stellar atmosphere code PHOENIX (Hauschildt & Baron 
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Table 8. 3?ca and ^caS interpolated linear regression 


(j a aMg 


^ca (Case B) 0.26 0.93 0.27 
3?ca5 (Case B) 0.35 1.19 0.30 

^ca (Case A) 0.23 0.48 0.48 
3ftca5 (Case A) 0.24 0.75 0.31 


and 3?ca5 linear regression on 
the observed snpernova spectra with 
sufficient wavelength coverage using 
the values interpolated to t = 0. 


Table 9. 3^5^ interpolated luminosity measure precision 




(Case B) 0.06 0.79 0.07 
(Case A) 0.04 0.35 0.12 


linear regression on the ob¬ 
served supernova spectra with suffi¬ 
cient wavelength coverage using the 
values interpolated to t = 0. 
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Fig. 11.— !fisi time evolution (top panel, Case B) and 3^5^ correlation with values interpo¬ 
lated to t = 0 (bottom panel). 
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Fig. 12.— 3^5*5 time evolution (top panel), '^sis corrected for the time dependence (bottom 
panel). 
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Fig. 13.— ^siss time evolution (top), 3 ^ 5^55 corrected for the time dependence (bottom). 
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Table 10. ^sis and 3^5*55 linear regression with t = 0 extrapolated values 


cr a aMg 


(Case B) 
'^siss (Case B) 
'^SiS (Case A) 
^siss (Case A) 


0.03 

0.67 

0.04 

0.04 

0.80 

0.05 

0.02 

0.49 

0.05 

0.04 

0.59 

0.06 


^'SlsiS and ^siss linear regression 
on observed supernovsewith sufficient 
wavelength coverage using the t = 0 
interpolated values 



Fig. 14.— ^ca for Case B. Black: Observed supernovae, Blue: day 20 W7 PHOENIX synthetic 
spectra, Pink: day 25 W7 PHOENIX synthetic spectra. 
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Fig. 15.— 3fJcaS for Case B. Black: Observed supernovse, Blue: day 20 W7 PHOENIX synthetic 
spectra, Pink: day 25 W7 PHOENIX synthetic spectra. 
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1999, 2004) of the parameterized SN la deflagration model W7 where we varied the total 
luminosity bolometric luminosity in the observer’s frame. We calculated one set at 20 days 
after explosion (the fiducial time of maximum in V) and another at 25 days after explosion. 
We have assumed LTE for computational expedience. Studies by Baron et al. (1996) and 
Wheeler et al. (1998) have shown that LTE is a reasonable approximation in SNe la. The 
Si IIA6355 feature is not reproduced by W7, as shown by Baron et al. (2005), who used 
extremely detailed NLTE. Only the general trends of the synthetic spectra can then be 
expected to be matched by the real data. The use of PHOENIX with W7 models is described 
in detail in Lentz et al. (2001). 

We display in Eigs. 14-15 the comparison between ^ca and calculated on observed 
supernovaeof Table 3 (all dates) and on PHOENIX spectra. These synthetic spectra cover 
a large range in Mb- Since the blue magnitude increases monotonically in the PHOENIX 
models with the bolometric magnitude, these figures indicate the existence of a spectroscopic 
sequence with the bolometric luminosity. The slope of the observed spectra is closer to that 
obtained from the models 25 d after explosion. The increase in 3?caand ^cas begins for 
lower blue magnitudes in the PHOENIX synthetic spectra than in the observed spectra. 

Eigs. 16-18 display the same synthetic and observed spectra for ^si, ^SiS, and ^siss- 
Again, the existence of the spectroscopic sequence is clear in both the observed and the 
synthetic spectra. Moreover, while the general trend is only qualitatively similar for ^ca 
and 3?ca5, 3^5* and to a lesser extent ^siS and S^si^^show that these ratio behave in the 
same way with respect to the blue magnitude evolution for synthetic as well as for real 
spectra. 


8. Use of line ratios in the SNAP context 

In order to estimate the luminosity measure accuracy of these spectral indicators in a 
cosmological context, we simulated SNAP exposures (SNAP Collaboration 2005) at differ¬ 
ent redshifts. We used them to estimate the accuracy on the determination of the mean 
luminosity of “Branch normal” supernovae, which is the value needed to devise the Hubble 
diagram and the determination of cosmological parameters. 



29 




Fig. 16.— 3^5^ for —19.16 (top) and —17.58 (bottom) SN 1986G blue magnitude. Black: 
real supernovse. Blue: day 20 W7 PHOENIX synthetic spectra. Pink: day 20 W7 PHOENIX 
synthetic spectra. For Day 10 there is a discontinuity in Sft^jaround Mb = —18.9, this is 
real and is caused by the change in the formation of the Si llA5972 feature as discussed in 
Bongard et al.(in preparation). The very steep dependence of the various 3^5* ratios as of 
function of Mb in the range —19.5 < Mb < —20 is remarkably well reproduced by PHOENIX 
at day 25. The data points at low luminosity suggest a smaller slope. This justihes the 
limited luminosity range for the evaluation of the slopes in §§ 3-6 
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Fig. 17.— ^sis for —19.16 (top) and —17.58 (bottom) SN 1986G blue magnitude. Black: 
real supernovse. Blue: day 20 W7 PHOENIX synthetic spectra. Pink: day 20 W7 PHOENIX 
synthetic spectra. 















RSiSS RSiSS 


31 




Fig. 18.— ^siss for —19.16 (top) and —17.58 (bottom) SN 1986G blue magnitude. Black: 
real supernovse. Blue: day 20 W7 PHOENIX synthetic spectra. Pink: day 20 W7 PHOENIX 
synthetic spectra. 
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8.1. SNAP simulator 

8.1.1. The SNAP spectrometer 

To simulate the SNAP spectrometer, we assumed 50% optical and 70% CCD efficiencies 
(Kim et al. 2004). We considered the mirror snrface to be tt (i.e. a 2m mirror diameter) 
and implemented the re-binning of the snpernovaespectra to SNAP binning in the observer 
rest frame. The red and blue channel resolution of SNAP we used are summarized in 
Table 11 (A. Ealet, private communication; Ealet et al. [2003]). 

The separation between the two spectrograph channels is not yet precisely defined, 
therefore we treated the spectra as if the red and the bine detector properties were the same. 
Since the edges of the spectral range are harder to calibrate, their actual location can impact 
the final results. 


8.1.2. The simulated supernovce 

In order to simnlate SNAP exposnres, we calibrate each snpernova to its absolnte bine 
magnitnde and then redshift it to the desired 2 :. 


JV^M = iV,(0)(l + 2)(^) (8) 

Eq. (8) shows how the nnmber of photons per nnit time and per nnit area in the 

redshifted B band in the supernova rest frame at redshift N.y{z) = is 

related to ^^(0) = the number of photons per unit of time and of surface in 

the B band at rest in the supernova rest frame for a supernova at 10 pc. 

The finx dilntion factor for a snpernova at 2 : is A- where do = 10 pc, bnt the reddening 
factor 1 + z included in dL gets cancelled when photons are counted in redshifted spectral 
bands. We also introduced electronic and statistical noise using the values given in Table 12. 

The electronic noise for each 2000s exposnre was calcnlated nsing 

’ ^electronic 

and assumed to be a Gaussian with zero mean. The Poisson noise of the signal was also 
included for each of the 2000 s exposures needed to complete the total exposure time. 


= 17 ^ 


t 


2 X 103 


( 9 ) 
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Table 11. SNAP blue and red channel resolution 


BLUE RED 

MA) MA) A (A) a(A) ^ 


4000 

295. 

7500 

95 

4500 

235. 

8000 

90 

5000 

185. 

8500 

85 

5500 

155. 

9000 

80 

6000 

130. 

9500 

78 

6500 

115. 

10000 

75 

7000 

105. 




10000. 

78. 

14000. 

86. 

10500. 

77. 

14500. 

90. 

11000. 

77. 

15000. 

93. 

11500. 

78. 

15500. 

97. 

12000. 

79. 

16000. 

100. 

12500. 

80. 

16500. 

105. 

13000. 

82. 

17000. 

no. 

13500. 

84. 




“^The resolution is taken from A. Ealet (private commu¬ 
nication). 


Table 12. Noise per pixel per 2000 s. 


readout noise/pixel a/ 2 ■ 5e 

slow drift contribution 7 ■ 1.2e 

Leakage cnrrent flnct. for 2000s 6e 

Total electronic noise per pixel: lOe 

Presence of a cosmic ray (20% of cases) 12e 


Total Electronic noise per AA (3 spatial pixels) 17e 
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In Table 13 we display the number of photons per second expected in the B Bessel 
hlter for a SN laof absolute blue magnitude Mb = —19.42, not accounting for the CCD 
and optical transmission. The exposure time we use have been assumed to follow a (1 + z)^ 
power law for 1 < z < 1.7 and an exposure time of 10 hours for = 1.7, chosen in order to 
maintain a constant signal to noise as 2 : increases (A. Ealet, private communication). 


8.2. Results at z = 1.5 

We display in Fig. 19 the number of photons per bin simulated obtained for the standard 
supernova SN 1992A at z = 1.5 and 2 : = 1.7, respectively. This quantity differs from the 
usual flux per A, which explains the flatter shape of the spectrum as there is a difference 
of 1/A between them. For 2 = 1.7, the wavelength coverage will be too small to allow 
the calculation of ^si and ^sis, although it would be accessible to JEDI (the proposed 
spectroscopic wavelength coverage for JEDI is 0.8-2 /rm). 

We simulated 100 exposures for every supernova listed in Table 14 and calculated the 
spectral indicators previously defined for each of them and found the mean and standard 
deviation (see Tables 15, 16, 18, and 19). From this we calculate aMg assuming that the 
slope of the linear regression is known perfectly. In order to concentrate on the accuracy of 
the spectroscopic indicators in SNAP, and since at redshifts of 2 ~ 1 time dilatation doubles 
the evolution time of the light curve, we assumed that the spectra will be taken almost at 
maximum light. We thus included in our sample only the supernova spectra that were the 
close to the light curve maximum. 


8 . 2 . 1 . 

We list in Table 15 the luminosities measured for the simulated supernovaewithout any 
selection cut on the simulated 3^5^ values for Case A. The RMS dispersion due to the noise 
is displayed as well as the mean ^si used to calculate the blue magnitude. In this first table 
we used the linear regression calculated with 3^5^ values extrapolated to t = 0 when possible. 
Table 16 lists the same values but for the linear regressions corresponding to Case B. The 
large discrepancies observed for SN 1991T and SN 1991bg between the inferred (using 3?s'j) 
luminosity and the actual value suggests that the linear regression is invalid at low and high 
magnitudes, we discuss this further below; however, no Arm conclusion can be drawn from 
the single cases in our current sample. 

In order to estimate the mean accuracy for JDEM we add the values of SMb listed in the 
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Table 13. SNAP simulator s ^ in Bessel B 


z 

Blue filter {N^ s ^ cm^) 

Exposure Time 

1.0 

4.32 

5980s 

1.5 

1.75 

22700s 

1.7 

1.3 

36000s 


“^The number of photons does not take into ac¬ 
count the factor 0.5 *0.7 from the CCD efficiency. 
It is calculated here for a Mb = —19.42 SN la. 


Table 14. Supernovaesimulated in the SNAP context. 


Supernova 

Date (wrt maximum) 

Mb 

SN 1991T 

0 

-20.01 

SN 1998aq 

0,1,2 

-19.56 

SN 1981B 

0 

-19.54 

SN 1998bu 

-2,-1 

-19.49 

SN 1994D 

-3 

-19.47 

SN 1996X 

-2,0,1 

-19.43 

SN 1989B 

-1 

-19.42 

SN 1992A 

-1 

-19.32 

SN 1986G 

-1 

-17.58 or -19.16 

SN 1991bg 

0 

-16.62 
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Fig. 19.— SNAP SN 1992A 2 ; = 1.5 simulated exposure (top), and 2 ; = 1.7 (bottom). Black 
curves are only the supernova, red curves include Poisson and electronic noise. 
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last column of the table in qnadratnre and divide by iV — 2. We lose two degrees of freedom 
by finding the mean and the slope with the same sample of supernovaein Tables 15-16 (and 
then we assume that the slope used is known perfectly); these values are listed in Table 24. 
Red, dim SN 1991bg-like SNe la will be strongly selected against due to Malmquist bias in 
deep cosmological surveys, therefore we shall never consider them in our analyses (in Case 
A we inclnde SN 1986G in the 91bg-like class). Similarly SN 1991T-likes wonld also not be 
used for precision cosmology, but we quote results both including and excluding 91T-likes. 
In Case B where blue magnitude of SN 1986G is Mb = —19.16 as is favored by Reindl et al. 
(2005). the mean accnracy on the blue magnitude becomes ~ 0.17(0.66) blue magnitude 
excluding (inclnding) 91T-likes. 

Table 24 summarizes the total error expected in SNAP. We have calculated the total 
error by adding the “intrinsic Ums” (obtained from the linear regression fit) in quadrature 
with the (Tms estimated above due to noise from onr simnlations. We note that both error 
estimates include errors due to the same systematics and thus our estimate of the total error 
using Method 1 is qnite conservative 

We also note that our estimate of “the intrinsic aMg” strongly dependent on the slope 
of the linear regression and thus our restriction to Branch normals in obtaining the slope is 
optimistic. Nevertheless, the results of the synthetic spectral calculations bolster onr view 
that the linearity is more prononnced in the Branch normal region. For using the 
t = 0 extrapolated linear regression instead of the one without correction increases the slope 
enough to give a significant effect in the less favored Case A, by ~ 0.1 blue magnitude. Thus, 
to take full advantage of this luminosity measnre, we need spectra as close to to maximum 
light as possible, which will be feasible at z = 1.5 since the elapsed time between —1 and 
+1 days around maximum is 2.5 x 2 = 5 days and since SNAP will try to obtain spectra 
near peak to maximize signal to noise, the proposed JEDI mission will easily fulfill this goal 
since they will obtain multiple spectra every 7 days. 

As shown in Table 17, SNAP is expected to find ~ 107 SNe laat 2 : = 1.5. The accuracy 
expected on the luminosity measure with 3^5^ would then be statistically improved by a 
factor of vl07 ~ 10 (although as discussed above this floor will truly be determined by 
the systematic error). The mean SNe lablue magnitude at z = 1.5 that is the quantity 
of cosmological interest, will be determined with an accnracy of ~ (0.2)/10 ~ 0.02blue 
magnitude in Case A, excluding 91T-likes. Clearly in Case A if 91T-likes are included the 
method is not useful. It seems likely that 91T-likes will be identifiable by JDEM, so the 
smaller error is more reasonable. In Case B where its blue magnitude is —19.16, the mean 
accnracy wonld be ~ (0.01 — 0.07), depending on whether 91T-likes are included. Thus 
appears to be a good candidate for a snpplementary tool to determine the eqnation of state 
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of the universe using supernovae. 


In Tables 18, and 19 we list the results of our simulations for both Cases A and B. 
Table 24 shows that ^siS is clearly a superior luminosity indicator in that the total error is 
nearly constant at ^ 0.2Mb for all cases. This makes this ratio appear to be an excellent 
secondary probe of SNe lablue magnitudes and to independently test evolutionary effects. 

In Case B, where 91T-likes are rejected, the preferred case, the mean blue magnitude 
accuracy becomes ~ 0.10 blue magnitude, including 91T-likes increases the dispersion to 
only 0.16 which is similar to the ~ 0.2 blue magnitudes dispersion of the SNe laafter stretch 
factor correction. We see that ^sis is more robust than ^si as the total error remains below 
0.20 blue magnitudes for any of our assumptions. 

We found that ^siss did not improve the results. This is due to the large wavelength 
binning in the SNAP design that already effectively transforms ^siS into an “integral” ratio 
compared to the lOA binning we previously used, as can be seen from the low dispersion due 
to Poisson noise(> 5%). Thus ^siss results are almost identical to those of ^siss, sometimes 
even a bit worse because the larger SNAP binning can include more of the Fe ll~ 5600A 
peak into the integral, reducing its quality as a luminosity indicator. 

With more supernovaeadded to our analysis it may be possible to include both 1991bg- 
like events into our calibration for the spectral indices. 


8.2.3. ^Ca 

Similar to ^siss, ^Cas has a smaller dispersion due to noise, but a larger calibration 
error because of the coarse SNAP binning. We thus only display ^ca results in this section. 
^Ca calibration also suffers from the coarse sampling of SNAP spectrograph, as could be 
seen in the rightmost column of Tables 20-21. Since we found no trend with epoch, we have 
used the results for the slope without correction to t = 0 (Table 4). 

As can be seen in these tables that summarize the luminosity measure precision esti¬ 
mated for 3?ca, the systematic error dominates, the error due to the Poisson noise being 
always lower than 0.19 blue magnitudes. 
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8.3. Evolution with 2 : 

SNAP wavelength coverage does not allow the nse either or ^sis for 2 : = 1.7. 
On the other hand 3?ca will be measnrable. JEDI wonld allow the measnrement of at 
z = 1.7, but 3?ca falls outside the range of its proposed IR spectrograph. Since the exposure 
time has been calculated so that the number of photons per bin will remain almost constant, 
the signal to noise will not be degraded too much. The dispersion due to noise (Tables 22-23 
and Table 24) remains small compared to the intrinsic dispersion of the supernovae. 

The ^ca blue magnitude accuracy remains dominated by the intrinsic dispersion of the 
supernovaecorresponding to (Tms ~ 0.2. This result can again be improved by using the 80 
supernovaeexpected. The blue magnitude precision would then be AMb ~ 0.02, which is 
within the SNAP cosmological requirements. 


9. Conclusions 

We have presented a way to automatically measure the value of the spectral indices ^ca 
and which should be robust. We have dehned three new spectral indices 3?cas, ^sis, 
and ^siss, which are less sensitive to noise in the spectrum. We have shown that for the 
limited set of Branch Normal supernovae where public spectra are available near maximum 
light that there is a good correlation between the spectral indices and the luminosity at peak. 
Within the caveat of small size of our sample we found that there is not a general trend in 
the spectral indices with epoch as long as the epoch is close to maximum light. We have used 
synthetic spectral models to show qualitatively that there does appear to be a deviation from 
linearity in the spectral indices, particularly among the dim, fast-declining supernovae, such 
as SN 1991bg. This result is in agreement with that of Phillips et al. (1999) who also found 
that there is a deviation from linearity in peak magnitude with 6 mi 5 as 617115 increases (see 
their Figure 8). It is also in agreement with some theoretical hydrodynamical results which 
suggest it may be difhcult to produce a SN la with such low nickel mass via a deflagration 
(W. Hillebrandt, private communication), although Hoflich et al. (2002) are able to vary 
the nickel mass by varying the density of the deflagration to detonation transition. It has 
been suggested that there is both a variation in SNe laluminosity with galaxy type (Branch 
et al. 1996; Hamuy et al. 2001) and that there are two populations of SNe laprogenitors 
(Scannapieco & Bildsten 2005; Mannucci et al. 2005) which could be related to the deviation 
in the linearity of the spectral indices. 

Restricting ourselves to the sample of Branch Normal supernovaewe hud that the pro¬ 
posed JDEM mission SNAP could in fact use ^sis as a complementary luminosity indicator. 



While the results may appear somewhat circular in that we have restricted our calibrators 
to have a small dispersion and then shown that the correlation based on that result also 
recovers the small dispersion, the ability to use a separate observable from light-curve shape 
will provide for a complementary measurement of the luminosity distance as a function of 
redshift in the proposed JDEM wide-held high-redshift spaced based snrveys snch as SNAP 
and JEDI. The lightcnrve based method of template htting will have a certain set of sys- 
tematics that depends on exactly how “like versus like” is chosen, whereas the spectral index 
method will have a different set of systematics. In addition, the spectral index methods are 
very weakly sensitive to both the effects of dnst and of the hnx calibration of the spectrnm, 
since they are dehned in a very narrow wavelength range. 

Our simulations show that the single spectrum obtained by SNAP for type identihcation 
will have sufficient signal to noise to measure the spectral indices to the required precision. 
JEDI, which will obtain mnltiple spectra for each snpernova, wonld be even better snited to 
nsing 'Stsis since they wonld have a larger sample of spectra very close to maximnm light. 

As more spectra from nearby supernovaewith known peak luminosities become available 
we should be able to refine our calibration, and thus improve its accuracy for JDEM and 
other snpernovaesnrveys. We note that onr actnal nnmerical resnlts rely on both a very small 
sample and also on the valne of the absolnte magnitndes assigned by Reindl et al. (2005), 
althongh since we qnote resnlts with both Case A and B the sensitivity to the assigned 
absolute magnitudes can be gauged somewhat. All available distance estimates (and hence 
absolute magnitudes) to our sample will contain errors, and we have chosen to use the 
results of Reindl et al. (2005) becanse they are derived in a specific somewhat homogeneons 
manner. Thns, we regard this work as a “proof of concept” and the specific nnmerical valnes 
as our best current estimate. Clearly when 300 Hubble-fiow supernovae are available from 
current ground-based surveys (i.e. in the next 3-5 years), this analysis should be recomputed. 
Once accurate values are in hand, comparing absolute magnitude (determined from light¬ 
cnrve shape) and '^siS determined at high redshift may provide a nsefnl measnre of the 
evolntionary effects in the high-redshift sample. We discnss this farther in a forthcoming 
publication where we study the formation of the Si llA5972 line. 

We thank Anne Ealet for information abont the characteristics of SNAP and Peter 
Nngent for helpfnl discnssions. This work was snpported in part by NASA grants NAG5-3505 
and NAG5-12127, NSP grants AST-0204771 and AST-0307323, and the Region Rhone-Alpes. 
PHH was supported in part by the Pole Scientifique de Modelisation Numerique at ENS- 
Lyon. This research nsed resonrces of: the National Energy Research Scientific Compnting 
Genter (NERSC), which is snpported by the Office of Science of the U.S. Department of 
Energy under Contract No. DE-AC03-76SE00098; and the Hochstleistungs Rechenzentrum 



Nord (HLRN). We thank both these institutions for a generous allocation of computer time. 



Table 15. 


3?5i (Case A) 


SN name 

mean 


Inferred Mb 

Mb 

6 Mb 

91T 

1.16 

0.81 

-16.84 

2.34 

-3.17 

98aq 

0.18 

0.08 

-19.66 

0.23 

0.10 

81B 

0.33 

0.09 

-19.24 

0.26 

-0.3 

98bu 

0.17 

0.05 

-19.70 

0.16 

0.21 

94D 

0.29 

0.09 

-19.35 

0.27 

-0.12 

96X 

0.27 

0.09 

-19.42 

0.27 

-0.01 

89B 

0.36 

0.09 

-19.14 

0.33 

-0.27 

92A 

0.26 

0.11 

-19.44 

0.36 

0.12 

86G 

0.32 

0.16 

-19.26 

0.47 

1.68 

91bg 

0.56 

0.37 

-18.58 

1.09 

1.96 


^6 Mb is the difference between the trne valne and 
that inferred nsing onr correlation with the spectral in¬ 
dex SfJs'j 










Table 16. 3?5i (Case B) 


SN name 

mean 


Inferred Mb 

O' Mb 

5Mb 

91T 

1.16 

0.81 

-18.33 

1.02 

-1.68 

98aq 

0.18 

0.08 

-19.56 

0.10 

0.00 

81B 

0.33 

0.09 

-19.37 

0.12 

-0.16 

98bu 

0.17 

0.05 

-19.58 

0.07 

0.09 

94D 

0.29 

0.09 

-19.43 

0.12 

-0.04 

96X 

0.27 

0.09 

-19.46 

0.12 

0.02 

89B 

0.36 

0.09 

-19.34 

0.15 

-0.08 

92A 

0.26 

0.11 

-19.47 

0.16 

0.14 

86G 

0.32 

0.16 

-19.39 

0.21 

0.22 

91bg 

0.56 

0.37 

-19.08 

0.48 

2.47 


^5Mb is the difference between the true value and 
that inferred using our correlation with the spectral in¬ 
dex ^si 


Table 17. Number of supernovaein redshift bins for SNAP. 















z 

0.6 

0.7 

0.8 

0.9 

1.0 

1.1 

1.2 

1.3 

1.4 

1.5 

1.6 

1.7 

# of SN 

150 

171 

183 

179 

170 

155 

142 

130 

119 

107 

94 

80 


^These results were obtained from (Aldering et al. 2004) 













Table 18. 3^5*5 Case A 


SN name 

mean 

<^5RsiS 

Inferred Mb 

O' Mg 

6 Mb 

91T 

0.66 

0.03 

-19.71 

0.07 

-0.30 

98aq 

0.69 

0.03 

-19.64 

0.06 

0.08 

81B 

0.78 

0.04 

-19.47 

0.07 

-0.07 

98bu 

0.86 

0.05 

-19.31 

0.11 

-0.18 

94D 

0.77 

0.04 

-19.48 

0.07 

0.01 

96X 

0.74 

0.04 

-19.54 

0.08 

0.11 

89B 

0.83 

0.06 

-19.37 

0.13 

-0.05 

92A 

0.89 

0.07 

-19.24 

0.15 

-0.08 

86G 

1.06 

0.07 

-18.89 

0.15 

1.31 

91bg 

0.94 

0.21 

-19.15 

0.43 

2.52 


^5Mb is the difference between the trne valne and 
that inferred using our correlation with the spectral 
index "^sis 










Table 19. ^sis Case B 


SN name 

mean 


Inferred Mb 

O'Mb 

6Mb 

91T 

0.66 

0.03 

-19.64 

0.05 

-0.36 

98aq 

0.69 

0.03 

-19.60 

0.04 

0.04 

81B 

0.78 

0.04 

-19.48 

0.05 

-0.06 

98bu 

0.86 

0.05 

-19.36 

0.08 

-0.13 

94D 

0.77 

0.04 

-19.49 

0.05 

0.01 

96X 

0.74 

0.04 

-19.53 

0.05 

0.10 

89B 

0.83 

0.06 

-19.40 

0.10 

-0.02 

92A 

0.89 

0.07 

-19.31 

0.11 

-0.01 

86G 

1.06 

0.07 

-19.04 

0.11 

-0.12 

91bg 

0.94 

0.21 

-19.23 

0.32 

2.61 


^5Mb is the difference between the trne valne and that 
inferred nsing onr correlation with the spectral index 

3 ^ 5.5 


Table 20. ^ca Case A 2 : = 1.5 


SN name 

mean 


Inferred Mb 

O'Mb 

5Mb 

91T 

0.88 

0.03 

-19.81 

0.02 

-0.20 

81B 

1.20 

0.08 

-19.53 

0.07 

-0.01 

98bu 

1.20 

0.06 

-19.53 

0.05 

0.04 

94D 

1.28 

0.08 

-19.46 

0.07 

-0.01 

96X 

1.42 

0.08 

-19.35 

0.07 

-0.08 

89B 

1.06 

0.08 

-19.65 

0.07 

0.23 

92A 

1.28 

0.10 

-19.47 

0.08 

0.15 

91bg 

2.43 

1.71 

-18.48 

1.46 

1.86 


^5Mb is the difference between the true value and that 
inferred using our correlation with the spectral index 
3^Ca- 















Table 21. 3?ca Case B z = 1.5 


SN name 

mean 


Inferred Mb 

O' Mb 

5Mb 

91T 

0.88 

0.03 

-19.65 

0.01 

-0.36 

81B 

1.20 

0.08 

-19.49 

0.04 

-0.04 

98bu 

1.20 

0.06 

-19.49 

0.03 

0.00 

94D 

1.28 

0.08 

-19.46 

0.04 

-0.01 

96X 

1.42 

0.08 

-19.39 

0.04 

-0.04 

89B 

1.06 

0.08 

-19.56 

0.04 

0.14 

92A 

1.28 

0.10 

-19.46 

0.05 

0.14 

91bg 

2.43 

1.71 

-18.92 

0.8 

2.29 


^6 Mb is the difference between the true value and that 
inferred using our correlation with the spectral index 

^Ca- 


Table 22. ^ca Case A 2 : = 1.7 


SN name 

mean 


Inferred Mb 

OMb 

5Mb 

91T 

0.80 

0.03 

-19.87 

0.02 

-0.13 

81B 

1.32 

0.07 

-19.43 

0.06 

-0.11 

98bu 

1.26 

0.14 

-19.48 

0.12 

-0.00 

94D 

1.25 

0.09 

-19.49 

0.07 

0.02 

96X 

1.3 

0.05 

-19.45 

0.04 

0.02 

89B 

1.04 

0.08 

-19.67 

0.07 

0.25 

92A 

1.32 

0.12 

-19.42 

0.1 

0.1 

91bg 

2.39 

1.41 

-18.51 

1.2 

1.89 


^5Mb is the difference between the true value and that 
inferred using our correlation with the spectral index 
3^Ca- 

















Table 23. 3?ca Case B z = 1.7 


SN name 

mean 


Inferred Mb 

O' Mb 

5Mb 

91T 

0.80 

0.03 

-19.68 

0.01 

-0.33 

81B 

1.32 

0.07 

-19.43 

0.03 

-0.10 

98bu 

1.26 

0.14 

-19.47 

0.06 

-0.02 

94D 

1.25 

0.09 

-19.47 

0.04 

-0.00 

96X 

1.3 

0.05 

-19.45 

0.02 

0.02 

89B 

1.04 

0.08 

-19.57 

0.04 

0.15 

92A 

1.32 

0.12 

-19.43 

0.06 

0.11 

91bg 

2.39 

1.41 

-18.93 

0.66 

2.31 


^6 Mb is the difference between the trne valne and that 
inferred nsing onr correlation with the spectral index 
^s^s. 
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Table 24. Summary of Spectral Indicator Precision for SNAP 


Index a Mb ^Mb ^Mb 

(intrinsic) (noise) (total) 


3^Ca (z = 1-5, no 91T, Case A) 
^Ca (z = 1-5, with 91T, Case A) 
^Ca {z = 1.5, no 91T, Case B) 
^Ca {z = 1.5, with 91T, Case B) 
‘^Ca {z = 1.7, no 91T, Case A) 
'^Ca {z = 1-7, with 91T, Case A) 
‘^Ca {z = 1.7, no 91T, Case B) 
^Ca {z = 1.7, with 91T, Case B) 
(no 91T, Case A) 

'Sisi (with 91T, Case A) 

'Sisi (no 91T, Case B) 

(with 91T, Case B) 

(no 91T, Case A) 

(with 91T, Case A) 

'Sisis (no 91T, Case B) 

^SiS (with 91T, Case B) 


0.16 

0.14 

0.21 

0.16 

0.16 

0.23 

0.10 

0.10 

0.14 

0.10 

0.19 

0.21 

0.16 

0.16 

0.23 

0.16 

0.15 

0.22 

0.10 

0.11 

0.15 

0.10 

0.18 

0.20 

0.12 

0.24 

0.27 

0.12 

1.31 

1.32 

0.10 

0.14 

0.17 

0.10 

0.65 

0.66 

0.05 

0.11 

0.12 

0.05 

0.16 

0.17 

0.04 

0.09 

0.10 

0.04 

0.16 

0.16 


'^The methods used to estimate the error is described in the text. 
The method is very conservative and donble connts the systematic 
error. Clearly ^sis withont 91T-likes is the preferred method. 
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